Time-oscillating Lyapunov modes and auto-correlation functions for quasi-one- 
dimensional systems 
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The time-dependent structure of the Lyapunov vectors corresponding to the steps of Lyapunov 
spectra and their basis set representation are discussed for a quasi-one-dimensional many-hard-disk 
systems. Time-oscillating behavior is observed in two types of Lyapunov modes, one associated with 
the time translational invariance and another with the spatial translational invariance, and their 
phase relation is specified. It is shown that the longest period of the Lyapunov modes is twice as long 
as the period of the longitudinal momentum auto-correlation function. A simple explanation for 
this relation is proposed. This result gives the first quantitative connection between the Lyapunov 
modes and an experimentally accessible quantity. 
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In recent years, statistical mechanics based on chaotic 
dynamics has drawn much attention as a scheme to 
bridge between a deterministic microscopic dynamics and 
a probabiHstic macroscopic theory. In a chaotic system, 
the difference of two nearby phase space trajectories, the 
so called Lyapunov vector, diverges rapidly in time, so 
that a part of the dynamics is unpredictable and a sta- 
tistical treatment is required. To justify this scheme, 
much effort has been devoted to discovering the Unk be- 
tween macroscopic quantities, such as transport coeffi- 
cients, and chaotic quantities, such as the Lyapunov vec- 
tor and the Lyapunov exponent (the exponential rate of 
divergence or contraction of the magnitude of the Lya- 
punov vector) Some works have concentrated upon 
finding chaotic properties of many-particle systems, and 
some progress has been made in the search for global 
properties, like the conjugate pairing rule for the Lya- 
unov spectrum (the ordered set of Lyapunov exponents) 



The stepwise structure of the Lyapunov spectrum is 
a typical chaotic feature of many-particle systems 0. 
These steps in the spectrum are accompanied by global 
wavelike structures in the corresponding Lyapunov vec- 
tors, or Lyapunov modes 0, S S 0- The significance 
of this phenomenon is that this structure appears in the 
vectors associated with the Lyapunov exponents that are 
closest to zero, therefore it is connected with the slow 
macroscopic behavior of the system. Originally, it was 
observed in systems with hard-core particle interactions, 
but very recently numerical evidence for Lyapunov modes 
in a system with soft-core interactions was reported J|. 
Some a naly tical approaches, such as random matrix the- 
ory [lO|, [llj , kinetic theory [IjI and periodic orbit theory 
|l3|. have been used in an attempt to understand this 
phenomenon. 

A clue to understanding the stepwise structure of the 
Lyapunov spectrum and the Lyapunov modes is in the 
behavior of the Lyapunov vectors associated with the 
zero-Lyapunov exponents. In this scenario we claim 



that the origin of the Lyapunov modes is the transla- 
tion invariances and conservation laws which dominate 
the global behavior the system. For a two-dimensional 
system of N particles in periodic boundary conditions 
the Lyapunov vectors can be written as {5x, 5y, Spx, Spy) 
where 5x, 6y, 5px and Spy are the iV-vectors containing 
the X and y spatial components, and the x and y momen- 
tum components of the Lyapunov vector. For this sys- 
tem there are six degenerate zero-Lyapunov exponents 
and the Lyapunov modes are linear combinations of one 
of the basis sets {(A,0,0,0), (0,A,0,0), {px,Py, 0,0)}, 
or {(0,0,A,0), (0,0,0, A), {0,0, Px,Py)}, where is the 
A^-dimensional null vector, A is an A'^-dimensional vec- 
tor with all components equal to 1/^/N, and p^ and Py 
are the A^-dimcnsional vectors containing the x and y- 
components of the momentum of each particle 0. No- 
tice that the second basis set is the conjugate of the 
first, where {—6p,5q) is the conjugate of {5q,6p). The 
first set of basis vectors are associated with spatial and 
time translational invariance, the second basis set with 
total momentum and energy conservation. Therefore, 
in the zero-Lyapunov modes we observe that the sets 



{5xf}, {5yf>}, {5p^:>}, {^71, {5xf>/px„5yf>/py,} 



or {Sp'^xj IPxji5p''y^' /pyj} can have constant components 
independent of the particle number index j Q. If the 
zero-Lyapunov exponents and their Lyapunov vectors 
can be regarded as the zero-th step then we would ex- 
pect to see X dependent Lyapunov modes for fc-vector 
kn — 2mT/Lx constructed from the basis vectors 

(0, sin(fc„xj), 0, 0), (0, cos(fc„x,), 0, 0), (1) 

(sin(fc„a;j), 0, 0, 0), {cos{Kxj), 0, 0, 0), (2) 
{pxj sm{knXj),Pyj sm{knXj),0, 0), 

{Pxj cos{knXj),Pyj cos{knXj), 0, 0), (3) 

in the quasi-one-dimensional system with periodic 
boundary conditions 6J. However, only the first two of 
these Q are observed (as the two-point steps or trans- 
verse modes) ^ ^J. The second two modes are longi- 
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tudinal and we notice that the last two basis vectors 
have time dependent normahsation coefficients. We 
can remove this time dependence by combining the time- 
translational invariance mode and the longitudinal modes 
in a particular way. It is these combined modes that form 
an approximate basis set for the numerically observed 
Lyapunov modes 

asm{ujnt){pxj sm{knXj),pyj sin(fc„a;j ), 0,0) 

+/3cos(w„t)(cos(fc„Xj), 0, 0, 0), (4) 

asin{uJnt){px3 COs{knXj),Pyj cos{knXj),0, 0) 

+/3 cos{uJnt){sin{knXj), 0, 0, 0), (5) 
acos{LUnt){pxj sm{knXj),Pyj sin(fc„a;j), 0, 0) 

+/3 sin(a;„t)(cos(fc„a;j), 0, 0, 0), (6) 

a cos{uint){px] COs{knXj),Pyj COs{knXj),0, 0) 

+(3 sm{ujnt)(sm{knXj),Q, 0, 0), (7) 

(where cun is a frequency) which are normalisable at all 
times. 

Although this scenario for the Lyapunov steps and 
modes seems to be plausible, convincing evidence is still 
absent. The wavelike structure in as a function of 

Xj, the so called the transverse spatial translational in- 
variance Lyapunov mode, is well known This mode is 
stationary in time, and categorizes one of the two types 
of Lyapunov steps. Ref. [fl] showed a mode structure in 
5y^J^^ /pyj as a function of x, namely the transverse time 
translational invariance Lyapunov mode. This mode ap- 
pears as a time-oscillating spatial wavelike structure, and 
categorizes another type of Lyapunov step. On the other 
hand, modes of the longitudinal components of the Lya- 
punov vectors are less convincing. Ref. claims a mov- 
ing mode structure in ^x^""* as a function of a;, in the 
same Lyapunov steps as the ones having the mode in 
^yf^ IPyjt but the relation between these two modes are 
not known. Another important unsolved problem is the 
time scale of the oscillation of the Lyapunov mode. The 
period is rather long, often thousands of times as long as 
the mean free time, and it should correspond to a macro- 
scopic or collective property of the system. However, no 
direct evidence of this has been reported. 

In this letter we show explicit numerical evidence for 
the time-oscillating wavelike structure of the multiple 
longitudinal Lyapunov modes, in particular Sx^"^ /p^j 

and ^x^""* as functions of Xj [basis vectors Q, (0), lO, 
(|7|]. and clarify the phase relation between them and 
the mode by'^^ jpyj. Further, we present numerical ev- 
idence to connect the time-oscillation of the Lyapunov 
modes to the time-oscillation of the momentum auto- 
correlation function. Oscillatory behavior has been ob- 
served for the momentum auto-correlation function in 
many macroscopic systems 0, Hsl. a nd has its origin in 
a collective motion of the systeinjlal • Our central result 
is that the longest period of oscillation of the Lyapunov 



modes is twice as long as the period of the longitudi- 
nal momentum auto-correlation function. Although the 
period changes with boundary conditions and with the 
number of particles A^, this inter-relation between peri- 
ods is always satisfied. A possible explanation for this 
relation is proposed. The auto-correlation function can 
be observed directly by the neutron and light scattering 
experiments [l^ I13 and is connected to transport coef- 
ficients using linear response theory so this result 
gives quantitative evidence of a connection between the 
Lyapunov modes and an experimentally accessible quan- 
tity. 

One of the difficulties of observing the Lyapunov steps 
and modes is that numerical calculations of Lyapunov 
spectra and vectors are very time-consuming, and analyt- 
ical methods for them are not well established, although 
there have been some attempts . On the other 

hand, if the above scenario for the Lyapunov steps and 
modes, based on universal properties like translational 
invariance and the conservation laws, can be justified, 
then a simple model should be sufficient to establish 
its origin. This is the motivation behind the quasi-one- 
dimensional many-hard-disk system proposed by 0, HO] ■ 
It is a very narrow rectangular system where the parti- 
cle ordering in the longer ai-direction is invariant. Parti- 
cle interactions are restricted to nearest neighbors only, 
so the numerical calculation is much faster than for the 
full two-dimensional system. Moreover it shows a much 
longer region of Lyapunov steps and clearer Lyapunov 
modes, compared with the full two-dimensional system. 
Another useful technique to get clear Lyapunov steps 
and modes, is to use hard- wall boundary conditions ja]. 
Although hard-wall boundary conditions break transla- 
tional invariance and momentum conservation, there is 
a simple relationship between the Lyapunov steps and 
modes for systems with different boundary conditions 
. In this letter we use the quasi-one-dimensional many- 
hard-disk system with hard-wall boundary conditions in 
the longer x direction and periodic boundary conditions 
in the shorter y direction [the boundary condition (H,P)]. 
This system exhibits both types of Lyapunov steps, and 
is only marginally slower numerically because only the 
two particles at the ends of the system can collide with 
the walls. We use units where the mass and disk radius 
R of each particle is 1 , the total energy is N and the di- 
mensions of the system are chosen as = l.bNLy + 2R 
and Ly = 2i?(l + 10^^), respectively. 

Figure ^ shows the stepwise structure of the normal- 
ized Lyapunov spectrum for the quasi-one-dimensional 
system consisting of 100 hard-disks with (H,P) bound- 
ary condition. In the inset to this figure we show the full 
normalized spectrum for the positive branch (the nega- 
tive branch is obtained by the conjugate pairing rule). In 
this system the spatial translational invariance and the 
total momentum conservation in the longitudinal direc- 
tion are violated, and the number of the zero-Lyapunov 



FIG. 1: Stepwise structure of the Lyapunov spectrum nor- 
malized by the largest Lyapunov exponent for the quasi- 
one- dimensional many-hard-disk system with (H,P) boundary 
condition. The number of particles is N = 100. The white- 
filled circles correspond to the stationary modes while the 
black filled circles are the time-oscillating modes. Inset: The 
positive branch of the normalized Lyapunov spectrum. 



exponents is 4. We recognize a clear stepwise structure 
in which two-point steps and one-point steps appear al- 
ternately. 

One-point steps in the Lyapunov spectrum in Fig. 
correspond to the stationary mode Sy^^^^ as a function of 
Xj [basis vectors (^]. This structure is well known as 
the transverse spatial translational invariant Lyapunov 
mode, so we omit further discussion here. Figure |5] 
shows that the Lyapunov modes corresponding to the 
two-point Lyapunov steps have time-oscillating wavelike 



mode structures in dx"-^^ /p. 



XJ, 



and 5x'f'^ as 



functions of Xj and the collision number nt with almost 
the same period Tiya ~ 16300. More details of this nu- 
merical calculation are given elsewhere 0, HH • By inves- 
tigating time-dependent Lyapunov modes for the other 
Lyapunov exponents [2ll | , it is shown that the spatial part 
of Lyapunov vector components 6x and 5y corresponding 
to the Lyapunov exponents in the n-th two-point step are 
approximately expressed using the basis vectors ^ and 
0. Notice however, that for (H,P) boundary conditions 
kn = mi/ Lx- Note that the time-translational invariance 
Lyapunov modes, that is the terms containing the mo- 
mentum in basis vectors (@J, (jSj, © and (O have the 
same structure in the same Lyapunov exponent, and are 
orthogonal to the longitudinal spatial translational in- 
variant Lyapunov components of the same basis vectors. 
Large fluctuations in the longitudinal time translational 
invariance Lyapunov mode, like in Fig. [2Ib), come from 
the terms for the longitudinal spatial translational in- 
variance Lyapunov modes. The spatial nodes of these 
Lyapunov modes are pinned at the hard- walls 0| . 

We now discuss the connection between the time- 
oscillations of Lyapunov modes and the momentum auto- 
correlation function. FigureOlshows the longitudinal mo- 
mentum auto-correlation normalized by its initial value 
as a function of the collision number nt in the same sys- 




FIG. 2: Local time averages of; 



{Sxf'^'^^ /p^j)t; and (c) (Syf^'^'' /pyj)t 



(a) (Sxy--"'),; (b) 

as functions of the col- 
lision number nt and the normalized position {xj)t/Lx of the 
x-component of the j-th particle. The corresponding Lya- 
punov exponent is indicated by the arrow in Fig. Q On the 
bottom of each graph is the projected contour plot of the 
three dimensional graph at the values —0.08 (dotted lines), 
(solid lines), and -f0.08 (broken lines). 



tem as that used in Figs. n]and[21 Here, an arithmetic 
average of the auto-correlation functions for 11 disks in 
the middle of the system (far from the hard walls) is 
taken. This auto-correlation function initially shows a 
simple exponential damping, as discussed elsewhere j2ll |. 
so that part is omitted in Fig. |3| A clear time-oscillating 
behavior is observed in Fig. 01 that is nicely fitted to 
a sinusoidal function times an exponential decay. This 
fitting gives the oscillating period Tacf ~ 8290, which 
suggests the approximate relation 



T, 



lya 



2T, 



acf 



(8) 



between the oscillating periods of the Lyapunov mode 
and the momentum auto-correlation function. Numeri- 
cally this relation is satisfied irrespective of the number 
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oscillating period of the longitudinal momentum auto- 
correlation function. The momentum auto-correlation 
can be measured through experiments like the neutron 
and light scattering experiments, so our result give a di- 
rect connection between the Lyapunov mode with an ex- 
perimentally accessible quantity. 
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FIG. 3: The time-oscillating part of the longitudinal momen- 
tum auto-correlation function Cx, normalized by its initial 
value, as a function of the collision number nt. The system 
is the same eis that used in Figs. and |5| Here, the line 
is a fit of the data to a sinusoidal function multiplied by an 
exponential function. 

of particles and boundary conditions . 

The relation (jSJ can be obtained from the following 
physical argument. Consider the Fourier space auto- 
correlation function < Px{t)* Px{0) > for the longi- 
tudinal component Pxit) of the momentum that oscil- 
lates with frequency LOacf, namely < Px{t)* Pxifi) >~ 
(j){t) exp{jwac/i}, where <j){t) is the damping factor and 
the suffix * indicates the complex conjugate. On the 
other hand we have a term for the longitudinal time 
translational invariance Lyapunov mode, like the first 
component of the basis vectors ((SJ) and iQ, as 6qx ^ 
'ijji{t)px{t) expliwiya^} with a frequency ojiya, where "01 (t) 
is an amplitude factor. We assume that if the quan- 
tity 5qx oscillates persistently in time, then its auto- 
correlation function < 5qx{t)* 5qx{Q) > should also os- 
cillate in time with the same frequency LUiya, namely 
< 6qx{t)* Sqx{0) >~ ip2{t) Gxp{iu!iyat} with a damp- 
ing factor ?/'2(i)- (Actually we can show that the auto- 
correlation function for the spatial longitudinal compo- 
nent of the Lyapunov vector oscillates with the same fre- 
quency as the Lyapunov vector component itself in the 
quasi-one-dimensional system jlj.) It follows from these 
equations and the assumption that ip2{t) e^cpliujiyat} ~ 
ipiit)* -ipiiO) (j){t) exp{i(tjQc/ - ^iva)t}, that ujiya ^ 
LOacf/'i-- The frequency is inversely proportional to the 
oscillation period, so we obtain relation 0. 

In conclusion, we have discussed time-oscillating Lya- 
punov modes and the momentum auto-correlation func- 
tion in quasi-one-dimensional many- hard-disk systems. 
We clarified the phase relation between components of 
multiple time-dependent Lyapunov modes, one coming 
from the longitudinal spatial translational invariance and 
another from the time translational invariance. We 
showed that the longest time-oscillating period of the 
Lyapunov modes is almost twice as long as the time- 
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